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Abstract. We reconsider the basic properties of ray-transfer matrices for first- 
order optical systems from a geometrical viewpoint. In the paraxial regime of 
scalar wave optics, there is a wide family of beams for which the action of a ray- 
transfer matrix can be fully represented as a bilinear transformation on the upper 
complex half-plane, which is the hyperbolic plane. Alternatively, this action can 
be also viewed in the unit disc. In both cases, we use a simple trace criterion that 
arranges all first-order systems in three classes with a clear geometrical meaning: 
they represent rotations, translations, or parallel displacements. We analyze in 
detail the relevant example of an optical resonator. 
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1. Introduction 

Matrix methods [1,2] offer the great advantage of simplifying the presentation of linear 
models and clarifying the common features and interconnections of distinct branches 
of physics [3] . Modern optics is not an exception and a wealth of input-output relations 
can be compactly expressed by a single matrix [4] . For example, the well-known 2x2 
ray-transfer matrix, which belongs to the realm of paraxial ray optics, predicts with 
almost perfect accuracy the behavior of a Gaussian beam. 

In this respect, we note that there is a wide family of beams (including Gaussian 
Schell- model fields, which have received particular attention [5-15]) for which a 
complex parameter can be defined such that, under the action of first-order systems, 
it is transformed according to the famous Kogelnik ABCD law [16-20]. This is the 
reason why they are so easy to handle. This simplicity, together with the practical 
importance that these beams have for laser systems, explain the abundant literature 
on this topic [21,22]. 

The algebraic basis for understanding the transformation properties of such beams 
is twofold: the ray-transfer matrix of any first-order system is an element of the group 
SL(2, R) [23] and the complex beam parameter changes according to a bilinear (or 
Mobius) transformation [24]. 

The nature of these results seems to call for a geometrical interpretation. The 
interaction between physics and geometry has a long and fruitful story, a unique 
example is Einstein theory of relativity. The goal of this paper is precisely to provide 
such a geometrical basis, which should be relevant to properly approach this subject. 
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The material of this paper is organized as follows. In section [3 we include a brief 
review of the transformation properties of Gaussian beams by first-order systems, 
introducing a complex parameter Q to describe the different states as points of 
the hyperbolic plane. The action of the system in terms of Q is then given by 
a bilinear transformation, which is characterized through the points that it leaves 
invariant. From this viewpoint the three basic isometries of this hyperbolic plane 
(i.e., transformations that preserve the distance), namely, rotations, translations, and 
parallel displacements, appear linked to the fact that the trace of the ray-transfer 
matrix has a magnitude lesser than, greater than, or equal to 2, respectively. 

In section [3] we present a mapping that transforms the hyperbolic plane into the 
unit disc (which is the Poincare model of the hyperbolic geometry) and we proceed 
to study the corresponding motions in this disc. Finally, as a direct application, 
in section 0] we treat the case of periodic systems, which are the basis for optical 
resonators, providing an alternative explanation of the standard stability condition. 

We emphasize that this geometrical scenario does not offer any advantage in 
terms of computational efficiency Apart from its undeniable beauty, its benefit lies 
in gaining insights into the qualitative behaviour of the beam evolution. 

2. First-order systems as transformations in the hyperbolic plane H 

We consider the paraxial propagation of light through axially symmetric systems, 
containing no tilted or misaligned elements. The reader interested in further details 
should consult the extensive work of Simon and Mukunda [25-29] . We take a Cartesian 
coordinate system whose Z axis is along the axis of the optical system and represent 
a ray at a plane z by the transverse position vector x{z) (which can be chosen in the 
meridional plane) and by the momentum p(z) = n(z)dx/dz [30]. Here n(z) is the 
refractive index and dx/dz is the direction of the ray through z. 

At the level of ray optics, a first-order system changes the ray parameters by the 
simple transformation [31] 



where the primed and unprimed variables refer to the output and input planes, 
respectively, and M is the ray-transfer matrix that must satisfy the condition [32] 



which means that M is an element of the group SL(2, R) of real unimodular 2x2 
matrices. 

When one goes to paraxial-wave optics, the beams are described in the Hilbert 
space L 2 of complex- valued square- integrable wave- amplitude functions ip(x). The 
classical phase-space variables x and p are now promoted to self-adjoint operators by 
the procedure of wavization [33] , which is quite similar to the quantization of position 
and momentum in quantum mechanics. 

We are interested in the action of a ray-transfer matrix on time-stationary fields. 
We can then focus the analysis on a fixed frequency u, which we shall omit henceforth. 
Moreover, to deal with partially coherent beams we specify the field not by its 
amplitude, but by its cross-spectral density. The latter is defined in terms of the 
former as 




(2.1) 




detM = AD — BC = 1, 



(2.2) 



T(x 



i,x 2 ) = (ip*(x 1 )ip(x 2 )), 



(2.3) 
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where the angular brackets denote ensemble averages. 

There is a wide family of beams, known as Schell-model fields, for which the 
cross-spectral density (|2.3[) factors in the form 

r(a;i,a: 2 ) = [I(x 1 )I{x 2 )] 1/2 fi(x 1 ~ x 2 ), (2.4) 

where I is the intensity distribution and fi is the normalized degree of coherence, which 
is translationally invariant. When these two fundamental quantities are Gaussians 



I(x) = ■ — exp i --tj 



(2.5) 



H(x) = exp ( -^j 



x 2 



the beam is said to be a Gaussian Schcll model (GSM). Here X is a constant 
independent of x that can be identified with the total irradiance. Clearly, 07 and <7 M 
are, respectively, the effective beam width and the transverse coherence length. Other 
well-known families of Gaussian fields are special cases of these GSM fields. When 
071 <C 07 we have the Gaussian quasihomogeneous field, and the coherent Gaussian 
field is obtained when — > 00. In any case, the crucial point for our purposes is the 
observation that for GSM fields one can define a complex parameter Q [34] 

Q = l^ + l T^—^ ( 2 - 6 ) 
R koi a 

where 

- = - + ^ (27) 

and R is the wave front curvature radius. This parameter fully characterizes the beam 
and satisfies the Kogelnik ABCD law; namely, after propagation through a first-order 
system, the parameter Q changes to Q' via 

«' = *l M -« = OT < 2 - 8 > 

Since Im Q > by the definition l|2.6l) , one immediately checks that Im Q' > and we 
can thus view the action of the first-order system as a bilinear transformation on 
the upper complex half- plane. When we use the metric ds = \dQ\/ ImQ to measure 
distances, what we get is the standard model of the hyperbolic plane H [35]. This 
plane H is invariant under bilinear transformations. 

We note that the whole real axis, which is the boundary of H, is also invariant 
under l|2.8() and represents wave fields with unlimited transverse irradiance (contrary 
to the notion of a beam). On the other hand, for the points in the imaginary axis we 
have an infinite wave front radius, which defines the corresponding beam waists. The 
origin represents a plane wave. 

Bilinear transformations constitute an important tool in many branches of 
physics. For example, in polarization optics they have been employed for a simple 
classification of polarizing devices by means of the concept of eigenpolarizations of the 
transfer function [36,37]. 

In our context, the equivalent concept can be stated as the beam configurations 
such that Q = Q' in equation l|2.8|l . whose solutions are 
1 



2B 



(D-A)±^/{A + D) 2 -A . (2.9) 
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These values of Q are known as the fixed points of the transformation. 

The trace of M, Tr(M) = A + D, provides a suitable tool for the classification 
of optical systems [38]. It has also played an important role in studying propagation 
in periodic media [39]. When [Tr(M)] 2 < 4 the action is said elliptic and there are 
no real roots: they are complex conjugates and only one of them lies in H, while the 
other lies outside. When [Tr(M)] 2 > 4 there are two real roots (i.e., in the boundary 
of H) and the action is hyperbolic. Finally, when [Tr(M)] 2 = 4 there is one (double) 
real solution and the system action is called parabolic. 

To proceed further let us note that by taking the conjugate of M with any matrix 
C e SL(2,R) 

Mc = CMC -1 , (2.10) 

we obtain another matrix of the same type, since Tr(M) = Tr(Mc). Conversely, if two 
systems have the same trace, one can always find a matrix C satisfying equation (|2.10() . 

Note that Q is a fixed point of M if and only if the image of Q by C (i.e., 
^[C, Q]) is a fixed point of Mc- In consequence, given any ray-transfer matrix M one 
can always find a C such that Mc takes one of the following canonical forms [40,41]: 

cos(tf/2) sin(tf/2) 
-sin(i?/2) cos(i?/2) 



Kc(0) = 



A c(0- ( „ e % ) - (2.11) 

TVT / N ( 1 

where < •& < 47r and £, v G M. These matrices define the one-parameter subgroups 
of SL(2, R) and have as fixed points +i (elliptic), and oo (hyperbolic), and oo 
(parabolic), respectively. They are the three basic blocks in terms of which any system 
action can be expressed. Clearly, Kc(i?) represents a rotation in phase space, Ac(£) 
is a magnifier that scales x up by the factor e^/ 2 and p down by the same factor, and 
Nc(y) represents the action of a thin lens of power v (i.e., focal length 1/v) [30]. 
For the canonical forms H2.11[) . the corresponding actions are 

, _ cos(i?/2)Q ~ sin(0/2) 
Q ~ sin(i?/2)Q + cos(t9/2)' 

Q' = e~*Q, (2.12) 



Q' = Q + v. 

The first is a rotation, in agreement with Euclidean geometry, since a rotation has 
only one invariant point. The second is a translation because it has no fixed points in 
H and the geodesic line joining the two fixed points (0 and oo) remains invariant (it 
is the axis of the translation). The third one is known as a parallel displacement. 

When one of the parameters 9, £, or v in Ij2.12|l varies, the transformed points 
Q' describe a curve called the orbit of Q under the action of the corresponding one- 
parameter subgroup. In figure l.a we have plotted typical orbits for the canonical forms 
H2.11f) . For matrices Kq(i?) the orbits are circumferences centered at the invariant 
point +i and passing through Q and —1/Q. For Ac(£), they are lines going from 
to the oo through Q and they are known as hypercicles. Finally, for matrices Nc(^) 
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_L (a) 
Elliptic Hyperbolic Parabolic 




Figure 1. Plot of typical orbits in the hyperbolic plane H: (a) canonical transfer 
matrices as given in equation 12.111 and (b) arbitrary transfer matrices obtained 
by matrix conjugation as in equation I2.1UI , 



the orbits are lines parallel to the real axis passing through Q and they are known as 
horocycles [42]. 

For a general matrix M the corresponding orbits can be obtained by transforming 
with the appropriate matrix C the orbits described before. The explicit construction 
of the family of matrices C is not difficult: it suffices to impose that C transforms 
the fixed points of M into the ones of Kc(tf), Ac(£), or Nc(f), respectively. Just to 
work out an example that will play a relevant role in the forthcoming, we consider a 
matrix M representing an elliptic action with one fixed point denoted by Qf . Since 
the fixed point for the corresponding canonical matrix Kc(t?) is +i, the matrix C wc 
are looking for is determined by 

9[C,Q f ]=i. (2.13) 

If the matrix C is written as 

C\ C2 



C3 C4 



(2.14) 



the solution of l|2.13|l is 



CiRcQ f + C 3 ImQ f 

\Q f \ 2 

(2.15) 

_ dhnQf-CaBeQf 
04 — — 



\Qf\ 2 

In addition, the condition det C = +1 imposes 



that, together (|2.15|l determines the matrix C in terms of the free parameter C 
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In figure l.b we have plotted typical examples of such orbits for elliptic, 
hyperbolic, and parabolic actions. We stress that once the fixed points of the ray- 
transfer matrix are known, one can ensure that Q' will lie in the orbit associated to 
Q. 



3. First-order systems as transformations in the Poincare unit disc P 

To complete the geometrical setting introduced in the previous Section, we explore 
now a remarkable transformation (introduced by Cayley) that maps bijectively the 
hyperbolic plane H onto the unit disc, denoted by D. This can be done via the 
unitary matrix 

i 



in such a way that 



At = U M U 1 = ( " f, ) , (3.2) 



(i* a* 

where At is a matrix with det At = +1 and whose elements are given in terms of 
those of M by 

a=±[(A + D)+i(C-B)}, 

(3.3) 

(3=±[(B + C)+i(D-A)]. 

In other words, the matrices At belong to the group SU(1, 1), which plays an essential 
role in a variety of branches in physics. Obviously, the bilinear action induced by these 
matrices is 

B*+a*0 

where Q is the point transformed by (|3.1|) of the original Q: 

Q = T~h- (3-5) 

The transformation by U establishes then a one-to-one map between the group 
SL(2, M) of matrices M and the group SU(1, 1) of complex matrices At, which allows 
for a direct translation of the properties from one to the other. 

It is easy to see that HI maps onto D, as desired. The imaginary axis in EI goes 
to the Y axis of the disc ID (in both cases, R — oo and define beam waists). In 
particular, Q = +i is mapped onto Q — 0. The boundary of H (the real axis) goes 
to the boundary of D (the unit circle) , and both boundaries represent fully unlimited 
irradiance distributions (i.e., non-beam solutions). 

Since the matrix conjugation l|3.2|) does not change the trace, the same geometrical 
classification in three basic actions still holds. In fact, by conjugating with U the 
canonical forms l|2.11|l . we get the corresponding ones for SU(1, 1): 

exp(ztf/2) 

exp(-M?/2) 



a i< \ - I cosh (£/2) i sinh(£/2) 

~ I -i s inh(C/2) cosh^/2) ' ' 




Figure 2. Plot of typical orbits in the Poincare unit disc D: (a) canonical transfer 
matrices as given in equation 13.61 and (b) arbitrary transfer matrices. 



„/2 l + iv/2 ) ' 

that have as fixed points (elliptic), +i and —i (hyperbolic) and +i (parabolic), 
respectively. The first matrix represent a rotation in phase space, also called 
a fractional Fourier transformation, while the second one is sometimes called a 
hyperbolic expander [32]. 

The corresponding orbits for these matrices are defined by 



As plotted in figure 2. a, for matrices /Ce(i?) the orbits are circumferences centered 
at the origin. For J[c(£), they are arcs of circumference going from the point +i to 
the point — i through Q. Finally, for the matrices Afc{v) the orbits are circumferences 
passing through the points i, Q, and — Q*. In figure 2.b we have plotted the 
corresponding orbits for arbitrary fixed points. 



4. Application to optical resonators 

The geometrical ideas presented before allows one to describe the evolution of a GSM 
beam by means of the associated orbits. As an application of the formalism, we 
consider the illustrative example of an optical cavity consisting of two spherical mirrors 
of radii R\ and i?2, separated a distance d. The ray-transfer matrix corresponding to 




Figure 3. Plot of | Tr(M)| in terms of the parameters gi and 32 of the optical 
resonator. The plane | Tr(M)| = 2 is also shown. The density plot of the three- 
dimensional figure appears at the bottom. 



a round trip can be routinely computed [4] 
M 



^ 2^132 - 9i + 92 ~ 1 7j(2ffi52 + .91 + 92) 



2 

(2^152 -9i- 92) ^9\92 + 5i - .92 - 1 



(4.1) 



V d 

where we have used the parameters (i — 1,2) 

d , \ 

9i = 1 - p". (4-2) 

Note that 

Tr(M) = 2{2g l92 - 1). (4.3) 

Since the trace determines the fixed point and the orbits of the system, the 9 
parameters establish uniquely the geometrical action of the resonator. To clarify 
further this point, in figure 3 we have plotted the value of |Tr(M)| in terms of g± 
and 52- The plane | Tr(M)| = 2, which determines the boundary between elliptic and 
hyperbolic action, is also shown. At the top of the figure, a density plot is presented, 
with the characteristic hyperbolic contours. 

Assume now that the light bounces N times through this system. The overall 
transfer matrix is then M N , so all the algebraic task reduces to finding a closed 
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Figure 4. Plot of the successive iterates for typical elliptic, hyperbolic, and 
parabolic actions. In (a) the points are plotted in the hyperbolic plane H, while 
in (b) they are represented in the unit disc D. For hyperbolic and parabolic 
actions, the iterates tend to the real axis and the unit circle, respectively. 



expression for the Nth. power of the matrix M. Although there are several elegant 
ways of computing this power [21], we shall instead apply our geometrical picture: the 
transformed beam is represented by the point 

Q N = V[M,Q N _ 1 ] = *[M N ,Q ], (4.4) 

where Qo denotes the initial point. 

Note that all the points Qn lie in the orbit associated to the initial point Q by 
the single round trip, which is determined by its fixed points: the character of these 
fixed points determine thus the behaviour of this periodic system. By varying the 
parameters g of the resonator we can choose to work in the elliptic, the hyperbolic, or 
the parabolic case [43]. 

To illustrate how this geometrical approach works in practice, in figure 4. a we have 
plotted the sequence of successive iterates obtained for different kind of ray-transfer 
matrices, according to our previous classification. In figure 4.b we have plotted the 
same sequence but in the unit disc, obtained via the unitary matrix U. 

In the elliptic case, it is clear that the points Qn revolve in the orbit centered at 
the fixed point and the system never reaches the real axis. Equivalently, the points 
Qjv never reach the unit circle. 

On the contrary, for the hyperbolic and parabolic cases the iterates converge to 
one of the fixed points on the real axis, although with different laws [44] . In the general 
context of scattering by periodic systems this corresponds to the band stop and band 
edges, respectively [45-49]. 

What we conclude from this analysis is that the iterates of hyperbolic and 
parabolic actions produce solutions fully unlimited, which are incompatible with our 
ideas of a beam. The only beam solutions are thus generated by elliptic actions and, 
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according with equation i|4.3|) , the stability criterion is 

0<|2 fl i5a-l| = |co8(i?/2)|<l, (4.5) 

where ~d is the parameter in the canonical form Kc in equation i|2.11[) . Such a condition 
is usually worked out in terms of algebraic arguments using ray-transfer matrices, 
although the final results apply exclusively to scalar wave fields. 

Finally, we stress that real cavities resonate with vector fields. The situation 
then is far more involved because the vector diffraction for (polarized) electric fields is 
more difficult to handle, even for systems with small Fresnel numbers and the ABCD 
law does not apply to the corresponding kernel [50] . Exact solutions for these vector 
beams have recently appeared [51]. In any case, there is abundant evidence that 
the stability condition 14.5fl works well. This could be expected, since the transition 
to scalar theories captures all the essential physics embodied in the more elaborated 
vector analogues [52]. 

5. Concluding remarks 

In this paper, we have provided a geometrical scenario to deal with first-order optical 
systems. More specifically, we have reduced the action of any system to a rotation, 
a translation or a parallel displacement, according to the magnitude of the trace of 
its ray-transfer matrix. These are the basic isometries of the hyperbolic plane HI and 
also of the Poincare unit disc D. We have also provided an approach for a qualitative 
examination of the stability condition of an optical resonator. 

We hope that this approach will complement the more standard algebraic 
techniques and together they will help to obtain a better physical and geometrical 
feeling for the properties of first-order optical systems. 
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